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£SJ . (1) We introduce random discrete Morse theory as a computational scheme to measure 

the complicatedness of a triangulation. The idea is to try to quantify the frequence of 
discrete Morse matchings with a certain number of critical cells. Our measure will depend 
on the topology of the space, but also on how nicely the space is triangulated. 
. (2) The scheme we propose looks for optimal discrete Morse functions with an elemen- 

tary random heuristic. Despite its naivete, this approach turns out to be very successful 
even in the case of huge inputs. 

(3) In our view the existing libraries of examples in computational topology are 'too 
easy' for testing algorithms based on discrete Morse theory. We propose a new library 
containing more complicated (and thus more meaningful) test examples. 
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^ ! 1 Introduction 

Libraries of objects for algorithm testing are extremely common in computational geometry. 
Their set-up requires particular care: If a library consists of objects 'too easy to understand', 
then basically any algorithm would score great on them, thus making it impossible for the 
researcher to appreciate the efficiency of her algorithms. Of course, agreeing on what examples 
^ \ should be regarded as 'easy' is a hard challenge, and how to quantify complicatedness is even 

harder. 

In the present paper, we focus on computational topology, which deals with simplicial 
complexes in an abstract manner, i.e., without prescribing a shape, a volume, or the dimen- 
sion of a Euclidean space in which they embed. We present a possible statistical approach, 
which we call random discrete Morse theory. The mathematical background relies on For- 
man's discrete Morse theory from 1998 [27\ 128] . which in turn builds on Whitehead's simple 
homotopy theory, developed around 1939 [73]. (Especially important is Whitehead's notion 
of collapsibility, which is a combinatorial strengthening of the contractibility property.) 

Our idea is to create a quantitative version of these two theories. For example, we would 
like to be able to tell not only if a complex is collapsible or not, but also 'how easy it is' 
to find a collapsing sequence. To give a mathematical basis to this intuition, we consider a 
random model where we perform elementary collapses completely at random. The likelihood 
to find a complete collapsing sequence this way, will measure how easy it is to collapse the 
complex. Although this likelihood is, in most cases, too difficult to compute, we can estimate 
it empirically in polynomial time. The following elementary heuristic takes also into account 
complexes that are not contractible. 
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Figure 1: The graph Aj. 
Algorithm: Random Discrete Morse 

Input: A d-dimensional (abstract, finite) simplicial complex C, given by its list of facets. 

(0) Initialize Co = c\ = ... = Cd = 0. 

(1) Is the complex empty? If yes, then STOP; otherwise, go to (2). 

(2) Are there free codimension-one faces? If yes, go to (3); if no, go to (4). 

(3) Elementary Collapse: Pick one free codimension-one face uniformly at random and 
delete it, together with the unique face that contains it. Go back to (1). 

(4) Critical Face: Pick one of the top-dimensional faces uniformly at random and delete it 
from the complex. If i is the dimension of the face just deleted, increment Cj by 1 unit. 
Go back to (1). 

Output: The resulting discrete Morse vector (co, c\, c%, . . . , Q). 

By construction, Cj counts the critical faces of dimension i. According to Forman |27j . 
any discrete Morse vector (co, c\, C2, . . . , Cd) is also the face vector of a cell complex homotopy 
equivalent to C. 

Definition 1 The discrete Morse spectrum of a (finite) simplicial complex C is the col- 
lection of all possible resulting discrete Morse vectors produced by the algorithm Random 
Discrete Morse together with the distribution of the respective probabilities. 

Example: Consider the graph Aj of Figure [T] above. As there are no free vertices in it, 
Random Discrete Morse picks an edge uniformly at random and deletes it. If the edge 
chosen is the central bridge (which happens with probability 4), the output discrete Morse 
vector is (2,3). If any other edge than the central one is chosen, the output vector is (1,2). 
The discrete Morse spectrum is therefore {|-(1,2), =-(2,3)}; or, shortly, {(1,2), (2,3)} if we 
simply want to list the vectors of the spectrum. 

The algorithm Random Discrete Morse requires no backtracking, and 'digests' the 
complex very rapidly. The output (1, 0, 0, . . . , 0) is a certificate of collapsibility. If the output 
is different from (1,0, ... ,0), the complex could still be collapsible with a different sequence 
of free-face deletions. Every optimal discrete Morse vector has some positive probability to 
appear as output of the algorithm; unfortunately, this probability can be arbitrarily small. In 
case optimality is not reached, the algorithm still outputs something meaningful, namely (as 
already mentioned) the /-vector of a cell complex homotopy equivalent to the given complex. 

Since the output arrives quickly, we can re-launch the program, say, 10000 times, possibly 
on separate computers (independently). The distribution of the obtained outcomes yields 
an approximation of the discrete Morse spectrum. By the so-called Morse inequalities, each 
output vector is componentwise larger or equal than the vector of Betti numbers (/?o, • • • , fid)- 
When the spectrum stays 'close' to the vector of Betti numbers, we can regard the triangu- 
lation to be easy. This allows an empirical analysis of how complicated the complex is. 
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Figure 2: The bipyramid. 

We point out that the problem of finding optimal discrete Morse functions (with as fewest 
critical cells as possible) is NP-hard [35j [Hj; even the decision problem of whether some given 
(connected, finite) simplicial complex is collapsible or not is NP-complete [69J. We therefore 
should not expect to immediately find optimal discrete Morse vectors for any input. Indeed, 
one can easily construct examples on which our (or similar) random heuristic performs poorly; 
see Appendix A. 

However, for many triangulations of huge size, our elementary random heuristic produces 
optimal discrete Morse functions in almost 100% of the runs of the program. This could 
be interesting in the future also for homology computations. Discrete Morse functions (for 
general cell complexes) are implicitly computed in several homology algorithms that are 
based on fast (co-)reduction techniques, like the packages CHomP [20J, RedHom [17], and 
Perseus [64J. 

The paper is structured as follows: First we give details of our algorithm (Section [2]) and 
compare it with previous approaches (Section [3|). Then we survey the existing topological 
and combinatorial lower bounds for optimal discrete Morse vectors (Section U]). Finally, we 
describe and examine a collection of examples coming from several different areas of topology 
(Section [5]). In our opinion, the resulting library (Appendix B) is a richer and more sensitive 
testing ground for implementations based on discrete Morse theory. 

2 Details of the algorithm and computational complexity 

In the following, we give a more explicit description of our random heuristic. 

The first thing we do is to build the Hasse diagram of the given simplicial complex C, 
which represents the incidence structure of the face poset of C; see Figure [3] for an example. 
It takes 0(d ■ I ■ T) steps to construct the Hasse diagram of a simplicial complex, in case 
the complex is given by its facets (or to be precise, by its vertex-facet incidences). Here d 
is the dimension of the input complex, T the total number of faces, and I the number of 
vertex-facet-incidences [37]; cf. also [38J. 

Once the upward Hasse diagram and the downward Hasse diagram are set up (see below), 
we deconstruct a copy of the upward Hasse diagram in every run of our program by deleting 
(randomly picked) critical faces or pairs of faces in case there are free faces. We illustrate 
this with a concrete example. 

Example 2 (The bipyramid) The 2-dimensional boundary of the bipyramid of Figure [2] 
has 6 triangles, 9 edges, and 5 vertices. We list the faces level-wise in lexicographic order 
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Figure 3: The Hasse diagram of the bipyramid with one critical triangle (234), one matching 
edge (34-134), and four free edges (13, 14, 23, 24) highlighted. 

and identify each face by a label k l denoting the fc-th face of dimension i in the respective 
lexicographic list. 



I 2 : [1,2,4], 2 2 : [1,2,5], 3 2 : [1,3,4], J, 2 : [1,3,5], 5 2 : [2,3,4], 6 2 : [2,3,5], 
I 1 : [1,2], 2 1 : [1,3], 3 1 : [1,4], 4 1 : [1,5], 5 1 : [2,3], 6 1 : [2,4], 7 1 : [2,5], 8 1 : [3,4], 9 1 : [3,5], 

1°: [1], 2°: [2], 3°: [3], [4], 5°: [5]. 



Next, we initialize the Hasse diagram. Hereby, the graph of the Hasse diagram is stored 
twice. In the upward Hasse diagram, we level-wise list all inclusions of i-dimensional faces in 
(i + l)-dimensional faces, 

i = l: 1^{1 2 ,2 2 }, 2 1 S{3 2 ,4 2 }, 3V{^,5 2 }, 4* / {2*,4*h 5 V {5 2 , 6% 
g ] /{^,5 2 }, 7 1 ^{2 2 ,6 2 }, S^iS 2 ^ 2 }, 9 1 ^{4 2 ,6 2 }, 

i = 0: l /^{l 1 ,2 1 ,3 1 ,4 1 }, 2° /■ {l\5\6\ 7 1 }, 3° /> {2\ 5\ 8 1 , 9 1 }, 
4°/>{3 1 ,6 1 ,8 1 }, 5° / ^{4 1 ,7 1 ,9 1 }, 

while in the downward Hasse diagram we level-wise list the (j — l)-dimensional faces that are 
contained in the j-dimensional faces, 

3=2: 2 2 \{1 1 ,4\7 1 }, 3 2 \ {2 1 , 3 1 , 8% 

4 2 \{2 1 ,4 1 ,9 1 }, 5 2 \{5^,^,^}, e^ls 1 ,? 1 ^ 1 }, 

j = V. l^U ^ }, 2 1 \{1°,3 }, 4 1 \U°,5°}, 5 1 \{2°,3 }, 

6 1 \{2°,4°}, 7 1 \{2° 1 5°} 1 8*\{S°, 4°}, 9 J \{3°, 5°}. 



Here, 3 1 /*■ {I 2 , 3 2 } is the short notation for the inclusion of the edge 3 1 : [1,4] in the two 
triangles I 2 : [1,2,4] and 3 2 : [1,3,4]. 

During each run, the downward Hasse diagram is maintained, while a copy of the upward 
Hasse diagram is updated after the removal of a critical face or of a pair consisting of a free 
face and the unique face it is contained in. The sequence of updating steps for the above 
example could be as follows: 
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0. compute downward and upward Hasse diagram 

1. initialize copy of the upward Hasse diagram 

2. free edges: none 

3. select random critical triangle: 5 2 : [2,3,4] 

4. update upward Hasse diagram: 

i = l: iV{^,n 2 1 /«{3 2 a 2 }, S 1 ^{1\3% 4^{2 2 a 2 }, 5 1 /> {6 2 }, 
d^U 2 }, 7 1 ^{2 2 ,6 2 }, 8 1 /{3 2 }, 9 1 ^{4 2 ,6 2 } 

5. free edges: 5 J ,d> J ,S J 

6. select random free edge: 8 1 : [3,4] paired with 3 2 : [1,3,4] 

7. update upward Hasse diagram: 

»=i: iv{i 2 ,n ^/{^}, svi^.n ^v{^ 2 ,^}, svn 

i = 0: l^il^S^S 1 ,^}, 2°S{1 1 ,5 1 ,6 1 ,7 1 }, 3° {2 1 , 5 1 , 9 1 }, 
4°f{3\6% 5V{^,7^'} 

8. free edges: 2 1 ,3 1 ,5 1 , 6 1 

9. . . . 

The downward Hasse diagram tells us precisely which parts of the upward Hasse diagram 
we have to update. For example, the choice of the critical triangle 5 2 : [2,3,4] forces us to 
update, via 5 2 \ {5 1 , 6 1 , 8 1 }, the inclusions of the edges 5 ,6 ,<§ J in the upward Hasse 
diagram (by removing the triangle 5 2 as including face). 

Triangulations of closed manifolds initially have no free faces. Thus, we start with an 
empty list of free faces and immediately remove a random critical face. For triangulations of 
manifolds with boundary or general simplicial complexes, we first have to initialize the list 
of free faces. (This extra effort in computation time can be seen by comparing the respective 
run times for the examples knot and nc_sphere in Table[3](see also Section [53]) : With 0.813 
seconds, the 3-ball knot takes slightly longer per round than the 3-sphere nc_sphere with 
0.470 seconds.) 

Whenever we are done with one level of the Hasse diagram, we initialize the set of free 
faces for the next level below. Besides updating the upward Hasse diagram in each round, 
we also keep track of 

• the current list of free faces (and update this list whenever we delete a critical face or 
a pair consisting of a free face and the unique face it is contained in), 

• the current discrete Morse vector (cq, c±, . . . , Cd) (which is initialized by (0, 0, ... , 0) and 
updated by incrementing Cj by one whenever a critical face of dimension % is selected). 

At the end of every round, the resulting discrete Morse vector (co, ci, . . . , Cd) is stored along 
with its number of appearances in the various rounds. Eventually, we output the list of all 
obtained discrete Morse vectors together with their frequencies. 

2.1 Implementation in GAP 

We implemented our random heuristic in GAP [30]. In particular, we used GAP operations 
on lists and sets to initialize and update Hasse diagrams and respective lists of free faces. 
Our implementation is basic and has roughly 150 lines of code. 

The largest complex (in terms of number of faces) we tested our program on has face vector 
/ = (5013,72300,290944,495912,383136,110880) [lj. For this triangulation of a 5-manifold 
different(!) from a 5-ball it took in total 60:17:33 + 21:41:31 h:min:sec to first build the Hasse 
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diagram and then to run the random heuristic once, respectively. As resulting discrete Morse 
vector we obtained (1,0,0,0,0,0); thus, this non-trivial 5-manifold is collapsible. 

We point out that there is considerable room for improvement with respect to computation 
time. First of all, the Hasse diagram of a complex can be stored in terms of (sparse) boundary 
matrices on which fast elimination steps represent elementary collapses; see Joswig [M] for a 
discussion. In addition, it is way faster to perform matrix operations in, say, C++ compared 
to elementary set operations in GAP. Nevertheless, if it comes to compute (respectively to 
simplify a presentation of) the fundamental group of a simplicial complex, then GAP provides 
efficient heuristics; cf. Section [5.101 

3 Comparison with other algorithms 

There are three main previous algorithmic approaches that aim to compute optimal dis- 
crete Morse functions for simplicial complexes, one by Joswig and Pfetsch [35j, one by Eng- 
strom [25] . and one by Lewiner, Lopes, and Tavares [44] (cf. also Lewiner [43] ). 

Tools that allow to improve discrete Morse functions were provided by Hersh [33]; and 
King, Knudson, and Mramor [39] discussed how to improve discrete Morse functions for 
geometric complexes in M 3 . 

3.1 The algorithm of Joswig and Pfetsch 

This deterministic algorithm, apart from a complete backtrack search, is currently the only 
available implementation that actually determines optimal discrete Morse functions for all 
inputs. In the Joswig-Pfetsch approach [35], the problem of finding an optimal discrete Morse 
function is translated into a maximal matching problem for the underlying graph of the Hasse 
diagram with an additional acyclicity condition |19[ [27] . The acyclic matching problem is 
then solved as an integer linear program. 

For various small instances, the Joswig-Pfetsch approach successfully produces optimal 
discrete Morse functions [35 j . A first case, however, for which the associated integer linear 
program was too large to handle is for the 16-vertex triangulation poincare [14] of the 
Poincare homology 3-sphere with / = (16, 106, 180,90). Joswig and Pfetsch interrupted the 
computation after one week (although they did not make use of the fact that at least six 
critical cells are necessary since the fundamental group of the Poincare homology 3-sphere 
is non-cyclic; cf. Section \5.9\i . For the same instance our heuristic found the optimal Morse 
vector (1,2,2, 1) within 0.02 seconds. Also for other small instances our heuristic was much 
faster, e.g., for Rudin's ball rudin [651 [74] with / = (14, 66, 94, 41), Joswig and Pfetsch needed 
103.78 seconds to achieve the optimal discrete Morse vector (1,0,0,0), while our heuristic 
found the optimum in 0.004+0.00107 seconds; cf. Section [531 

3.2 The approach by Engstrom 

The elegant approach by Engstrom |25] is deterministic and fast. Roughly speaking, the 
idea is to proceed by deleting vertex stars, rather than by deleting pairs of faces. Engstrom 
introduces what he calls 'Fourier-Morse theory', a theory based on Kahn-Saks-Sturtevant's 
notion of non-evasiveness, much like Forman's discrete Morse theory was based on White- 
head's notion of collapsibility. Instead of computing discrete Morse functions, Engstrom's 
algorithm computes optimal Fourier-Morse functions, which are some discrete Morse func- 
tions, but not necessarily the optimal ones among them. In particular, obtaining an output 
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(1,0,... ,0) with this algorithm yields a certificate of non-evasiveness, a stronger property 
than collapsibility; cf. [55] . 

However, there is a 3-ball with only 12 vertices which has the collapsibility property, but 
not the non-evasiveness property [12J. As for other examples, Engstrom obtains (1,5,5,1) 
as the optimal discrete Fourier-Morse vector for the 16-vertex triangulation of the Poincare 
homology 3-sphere poincare. Instead, the optimal discrete Morse vector for this example 
is (1,2,2,1). For Rudin's ball rudin Engstrom found (1,2,2,0) compared to the optimum 
(1, 0, 0, 0). So using this algorithm, the optimum is often missed, even on triangulations with 
relatively small size. 

3.3 The heuristic of Lewiner, Lopes, and Tavares 

The heuristic approach of Lewiner, Lopes, and Tavares [2] (cf. also Lewiner [33]) is fast and 
was used to produce optimal discrete Morse vectors for several large 2- and 3-dimensional 
complexes. The problem of finding optimal discrete Morse vectors is reformulated in terms 
of finding maximal hyperforests of hypergraphs. Then different greedy heuristics are used to 
obtain large hyperforests. 

It has to be remarked, though, that most of the instances listed in [H] and later in [33] 
are mostly harmless from the point of view of discrete Morse theory; they mainly are 2- 
dimensional surfaces or shellable 3-dimensional balls and spheres, or products thereof - 
with the exception of the three more complicated examples knot, nc_sphere, and bing from 
Hachimori's simplicial complex library [32]. It is precisely on these three examples that the 
greedy heuristics of Lewiner et al. produce somewhat inconsistent results. 

In [44], (1,1,1,0) was obtained for bing and knot. In |43j on p. 92, bing and knot 
appear with (1,0,0,0) without mentioning the improvement with respect to [44] . Moreover, 
nc_sphere is listed on p. 92 of |43] with (1,2,2,1) and it is noted on p. 89: "Trickier, the 
non-shellable 3-sphere (NC Sphere) is a delicate model since no discrete Morse function can 
reach the minimal number of critical points for smooth homotopy." This latter statement 
is false as we found (in 12 out of 10000 runs) the optimal discrete Morse vector (1,0,0,1) 
for nc_sphere; cf. Table [3l In fact, the 3-sphere nc_sphere with / = (381,2309,3856, 1928) 
is obtained from the 3-ball knot with / = (380,1929,2722,1172) by adding a cone over 
the boundary of knot. By this, every discrete Morse function on knot with discrete Morse 
vector (1, 0, 0, 0) can be used to produce a discrete Morse function with discrete Morse vector 
(1, 0, 0, 1) on nc_sphere. In contrast, it would theoretically be possible to have knot with 
optimal discrete Morse vector (1, 1, 1,0), while nc_sphere has optimal discrete Morse vector 
(1,0,0, 1). The best discrete Morse vector we found in 1000000 runs for knot is (1, 1, 1,0); 
see Table [3] — whereas, as mentioned above, Lewiner [43] seemed to claim (1,0,0,0) for this 
example, which would beat our algorithm. 

4 Theoretical lower bounds for discrete Morse vectors 

In this section, we briefly recall some theoretical lower bounds for minimal discrete Morse 
vectors. The obstructions for the existence of discrete Morse functions with a certain number 
of critical cells are of various nature. Here we basically use four different criteria. The 
first concerns ridge-facet incidences, the second follows from elementary algebraic topology 
(applied to the Morse complex), the third uses knot theory, and the fourth comes from smooth 
Morse theory. 
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4.0.1 Ridge- facet incidences and Euler characteristic 

In order for a collapse to start, there need to be free faces. This is how to create a first 
obstruction, namely by constructing (i-dimensional triangulations in which every (d — l)-face 
is contained in two or more (i-faces. 

The most famous example of this type is the dunce hat, a contractible 2-complex obtained 
from a single triangle by identifying all three boundary edges in a non-coherent way. In any 
triangulation of the dunce hat each edge belongs to either two or three triangles; cf. |llj . 
Hence, the dunce hat cannot be collapsible or, in other words, it cannot have (1,0,0) as 
discrete Morse vector. 

The vectors (1,0, 1) and (1, 1,0) are also forbidden for the dunce hat. In fact, since each 
elementary collapse deletes two faces of consecutive dimension, it does not change the Euler 
characteristic. In particular, the alternating sum of the entries of a discrete Morse vector 
should always be equal to the Euler characteristic of a complex. 

The dunce hat does, however, admit (1, 1, 1) as discrete Morse vector, which is therefore 
optimal. 

4.0.2 The Morse complex 

Forman showed that any discrete Morse vector on a simplicial complex C is also the face- 
vector of a model for C, that is, a CW-complex homotopy equivalent to C . 

Theorem 3 (Forman |28j) Assume that some d-complex C admits a discrete Morse func- 
tion with o L critical faces of dimension i (i = 0, . . . , d). Then C has a model with a i-cells, 
called Morse complex. 

This theorem results in several obstructions. First of all, the i-th (rational) Betti number 
of an arbitrary CW-complex is always bounded above by its number of i-dimensional cells. 

Corollary 4 (Forman's weak Morse inequalities |28j) Assume that some d-complex C 
admits a discrete Morse function with Ci critical faces of dimension i (i = 0, . . . , d). Then 
<k > Pi{C) for each i. 

The previous result still holds if we consider homology over a finite field. 

Corollary 5 Assume some d-complex C admits a discrete Morse function with Ci critical 
faces of dimension i (i = 0, . . . , d). Then Ci > dim Hi(C; Z p ) for each i and for each prime p. 

Sometimes it is convenient to focus on homotopy groups rather than on homology groups. 
Recall that the fundamental group of a CW-complex with one 0-cell is completely determined 
by its 2-skeleton; a presentation of the group can be obtained using the 1-cells as generators 
and the 2-cells as relators. In particular, if the CW-complex has no 1-cells, its fundamental 
group must be trivial; and if the CW-complex has only one 1-cell, its fundamental group 
must be trivial or cyclic. 

Corollary 6 Assume some d-complex C with fundamental group G admits a discrete Morse 
function with 1 critical face of dimension and c\ critical faces of dimension 1. Then 
c\ > rank(G), the minimal number of generators in a presentation of G. (In particular, if G 
is non-abelian, then c\ > 2.) 
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4.0.3 Knot-theoretic obstructions 

Obstructions coming from short knots have been considered first by Bing [T3], Goodrick |31j . 
and Lickorish [46], and later investigated by the two authors [91 1121 f52] and others. Recall 
that a knot K inside a triangulation of a 3-sphere is just a 1-dimensional subcomplex homeo- 
morphic to a 1-sphere (or in other words, a closed path in the 1-skeleton.) The knot group is 
the fundamental group of the knot complement inside the sphere. Knot groups are perhaps 
the main invariant in knot theory. 

In the simplest form (that is, for 3-dimensional spheres) the obstructions are of the fol- 
lowing type: 

Theorem 7 (Lickorish |46j; cf. also [9]) Assume some triangulated 3-sphere S admits 
some discrete Morse function with C2 critical 2-faces. Then, for any knot K inside S, one 
has 

c 2 > Tunk{G K ) - fi(K), 
where Gk is the knot group of K and fi(K) is the number of edges of K . 

The previous theorem is usually applied together with the following two well-known facts: 

(1) there are knots whose groups have arbitrarily high rank; for example, the knot group of 
a connected sum of m trefoils has rank > m + 1 (Goodrick [31j): 

(2) any knot can be realized with only 3 edges in a suitably triangulated 3-sphere (Bing |13j). 
In particular, if we consider a 3-sphere S containing the connected sum of three trefoils re- 
alized on three edges, then Theorem [7] yields C2 > 1 for all discrete Morse vectors (1, c±, C2, 1). 
Note that c\ = c^, because of Euler characteristic reasons. 

A similar statement can be proven for 3-dimensional balls. 

Theorem 8 ([9, Corollary 4.25]) Assume some triangulated 3-ball B admits some dis- 
crete Morse function with c\ critical edges. Let K be a knot in the 1-skeleton of B, realized 
as a path of b edges in the boundary of B plus a path of e = fi(K) — b interior edges. Then 

c\ > rank(Gx) — 2e, 

where Gk is the knot group of K. 

4.0.4 Morse-theoretical obstructions 

Very recently, the first author proved the following result for smooth manifolds. 

Theorem 9 (} 10\j ) Every smooth Morse vector is also a discrete Morse vector on some (com- 
patible) PL triangulation. In dimensions up to 7, the converse holds too. 

The converse statement is interesting for us because it yields further obstructions. For 
example, we know from the work by Boileau and Zieschang [15] and others, that for every 
r > 0, there is a (smooth) 3-manifold M r of Heegaard genus g > rank(Af r ) + r. It follows 
that for every PL triangulation T of M r , every discrete Morse vector on T has c\ > g > 
rank(M ra ) + r critical edges. 
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5 Towards a new library of triangulations 



Table provides a library of 45 instances for which we sampled the discrete Morse spectrum. 
We ran our random algorithm 10000 rounds on each example, except for eight examples for 
which we did fewer runs. The 45 examples were selected for different reasons as we will 
explain below. The respective examples are listed at the beginning of each subsection. The 
library of examples can be found online at |56| . 

An additional infinite series of complicated triangulations, based on a handlebody con- 
struction of Akbulut and Kir by [3] , was recently given in [70] . 

5.1 'Trivial' triangulations 

Examples: dunceJiat, d2nl2g6, regular_2_21_23_l 

Discrete Morse theory is trivial on 1-dimensional complexes (graphs) and 2-dimensional com- 
pact manifolds (surfaces); cf. [H]. A simple modification of our heuristic allows to incorporate 
this as follows. Once we reduced a simplicial complex to a 1-dimensional complex, we switch 
to a deterministic strategy: As long as there are edges that are contained in a cycle, delete 
such (critical) edges iteratively; then collapse the remaining tree/forest to a point/a collection 
of points. 

Definition 10 Let C be a connected finite simplicial complex. The normalized discrete 
Morse spectrum of C is obtained from the discrete Morse spectrum of C by normalizing 
every discrete Morse vector (co, c%, C2, . . . , Cd) in the spectrum to (1, c\ — Co + 1, C2, . . . , Cd) and 
adding up the probabilities for the original vectors that have the same normalization. 

Example: The graph Aj of Figure [1] has normalized discrete Morse spectrum {1-(1, 2) } 
or, for short, {(1, 2)}. 

Lemma 11 Every connected simplicial 1-complex K with n vertices and m > n — 1 edges 
has normalized discrete Morse spectrum {(l,m — n + 1)}. 

The homology vector in this case is H*(K) = (Z,Z m - n+1 ), so the weak discrete Morse 
inequalities (see Section |5T9|) are sharp. 

Lemma 12 Every triangulation K of a closed (connected) surface of Euler characteristics x 
has normalized discrete Morse spectrum {(1,2 — x, 1)}. More generally, the same holds for 
every strongly connected 2-complex K without free edges. 

Proof: Triangulations of surfaces are strongly connected. Hence, after the removal of 
only one critical triangle the remaining complex collapses onto a 1-dimensional complex, i.e., 
C2 = 1. The conclusion follows from the previous lemma and extends to all strongly connected 
2-complex K without free edges. □ 

The example d2nl2g6 in Table [3] is the unique vertex-transitive, vertex-minimal neighborly 
triangulation of the orientable surface of genus 6 [5], the example regular_2_21_23_l is a 
regular triangulation of the orientable surface of genus 15 with 21 vertices |49l Ch. 5]. 




Something can also be said on complexes with few vertices: 

Theorem 13 (Bagchi and Datta |7J) Every ^-acyclic simplicial complex with at most 7 ver- 
tices is collapsible. 
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Table 1: Total time to find optimal discrete Morse functions (in a single run) for each of the 
combinatorial 3-manifolds with up to 10 vertices. 



Vertices\Types 


s 3 


S^S 1 


S^S 1 


All 


Total time 
(in Min:Sec.Frac) 


5 


1 






1 


0.008 


6 


2 






2 


0.008 


7 


5 






5 


0.012 


8 


39 






39 


0.060 


9 


1296 


1 




1297 


3.836 


10 


247882 


615 


518 
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Corollary 14 Every ^-acyclic simplicial complex K with at most 7 vertices is extendably 
collapsible and therefore has trivial discrete Morse spectrum {(1,0,0)}. 

The 7-vertex bound is sharp; the triangulation duncejiat of the dunce hat is an 8-vertex 
example of a non-collapsible contractible complex. 

5.2 3-manifolds with up to ten vertices 

For all 250 359 examples in the catalog [53 J of triangulations of 3-manifolds with up to 10 ver- 
tices, optimal discrete Morse vectors were found by a single run of our program each; see 
TableEJ 

Theorem 15 All 250 359 examples of triangulated 3-manifolds with up to 10 vertices admit 
a perfect discrete Morse function. 

The spheres in this list are all shellable, as are all 3-spheres with up to 11 vertices |68j. 
The smallest known non-shellable 3-sphere S_3_13_56 (trefoil) has 13 vertices [52]. For all 
the 1134 non-spherical examples the statement of the theorem is new. 

5.3 Polytopal spheres 

Examples: S_3_100_4850, 600_cell, S_3_1000_2990 , S_5_100_472. 

We ran our program on the 3-dimensional boundary S_3_100_4850 of the cyclic 4-polytope 
with 100 vertices and 4850 facets, on the 3-dimensional boundary 600_cell of the 600-cell, 
on the 3-dimensional boundary S_3_1000_2990 of a stacked 4-polytope with 1000 vertices 
and 2990 facets, and on the 4-dimensional boundary S_5_100_472 of a stacked 5-polytope 
with 100 vertices and 472 facets. In all these cases we obtained the optimal discrete Morse 
vector (1,0,..., 0, 1) in 10000 out of 10000 tries. We also tested various other examples of 
simplicial polytopal spheres and we always observed a trivial spectrum in these experiments. 
However, the normalized discrete Morse spectrum of simplicial polytopal spheres is not trivial 
in general. 

Theorem 16 (Crowley, Ebin, Kahn, Reyfman, White, and Xue [22]) The 7 -simplex A 7 with 
8 vertices contains in its 2-skeleton an 8-vertex triangulation of the dunce hat onto which it 
collapses. 
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As a direct consequence of Theorem [TBI the 7-simplex A7 is not extendably collapsible. 
Therefore the spectrum of its boundary is non-trivial. The same can be said for the simplicial 
polytopal Griinbaum-Sreedharan 3-sphere No. 32, which contains the dunce hat [llj. 

5.4 Random spheres 
Example: S_3_50_1033. 

While random surfaces can easily be generated, we lack good random models for 3- or higher- 
dimensional manifolds; cf. [23]. One possible approach is to consider all triangulations of 
3-spheres or 3-manifolds with a fixed number n of vertices, with the uniform distribution. 
While this setting is very promising for performing random experiments, we need to get a 
hold on the set of all the triangulations with n vertices in the first place, a task that so far 
has only been solved for 3-manifolds with up to 11 vertices |68| . 

Another model can be derived by performing random walks on the set of all triangula- 
tions where each step is represented by a single bistellar flip. According to a theorem of 
Pachner [63], two distinct triangulations of a manifold are PL homeomorphic if and only 
if they can be connected by a sequence of bistellar flips. An implementation of bistellar 
flips for exploring the space of triangulations within one PL component is the program BIS- 
TELLAR [57]; see |14] for a program description. The bistellar flip approach for generating 
random triangulations depends on the number of executed flips as well as on the way the 
flips are chosen. As a consequence, triangulations with n vertices are not selected according 
to the uniform distribution. 

For the example S_3_50_1033, we started with the boundary of the cyclic 4-polytope with 
50 vertices and face vector / = (50, 1225,2350, 1175). We then applied 1500 bistellar 1-flips 
and reverse- 1-flips that were chosen randomly from all admissible flips. The resulting sphere 
S_3_50_1033 has /-vector (50, 1083, 2066, 1033). The average number of critical cells in 10000 
runs turned out experimentally to be roughly 3.2 (which is considerably larger than 2). We 
therefore can conclude heuristically that random spheres tend to have a non-trivial spectrum. 

5.5 Knotted triangulations of balls and spheres 

Examples: dunce_hat_in_3_ball, Barnette_sphere, B_3_9_18, trefoil_arc, trefoil, 
rudin, doublejtref oil_arc, doublejtref oil, triple_tref oil_arc, triple_tref oil, 
non_4_2_colorable, knot, nc_sphere, bing. 

The example dunce_hat_in_3_ball is a triangulated 3-ball that contains the 8-vertex trian- 
gulation dunceJiat in its 2-skeleton. To indeed get stuck with dunceJiat, we need to perform 
collapses without removing any of the 17 triangles of the dunce hat. This results in a low 
probability to get stuck. Indeed, in 1000 000 runs we always found (1,0,0,0) as resulting 
discrete Morse vector. 

The non-polytopal Barnette_sphere [8j with 8 vertices also has trivial observed spectrum: 
In 1000 000 runs of our program we obtained the optimal discrete Morse vector (1,0,0,1). 
For the non-shellable 3-ball B_3_9_18 [51J with 9 vertices and Rudin's non-shellable 3-ball 
rudin |65[ I74j with 14 vertices we achieved the optimal discrete Morse vector (1,0,0,0) in 
every run. Therefore, non-polytopality and non-shellability not necessarily cause a non-trivial 
observed spectrum. 

If we wish to construct triangulated balls or spheres of small size with a very non-trivial 
observed spectrum, we need to build in complicated substructures of small size (like compli- 
cated knots on few edges) to get stuck at. 
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Table 2: Average number of critical cells for the knotted spheres and their barycentric sub- 
divisions, based on 10000 random runs. 



trefoil 2.0778 double.tref oil 3.5338 triple_tref oil 5.9898 

trefoil_bsd 2.0202 double.tref oil_bsd 3.3414 triple.tref oil_bsd 5.7352 



The triangulated 3-sphere trefoil (S_3_13_56 |52j) contains a 3-edge trefoil knot in its 1- 
skeleton and has optimal discrete Morse vector (1, 0, 0, 1). This vector was obtained in roughly 
96% of the runs of our heuristic. The triangulated 3-sphere double_tref oil (S_3_16_92 
[12]) with optimal discrete Morse vector (1,0,0,1) has a 3-edge double trefoil knot in its 
1-skeleton. Here, (1, 0, 0, 1) was achieved only in 40% of the runs. The triangulated 3-sphere 
triple_tref oil (S_3_18_125 |12| ) contains a 3-edge triple trefoil knot in its 1-skeleton and 
has optimal discrete Morse vector (1, 1, 1, 1), which we found 30% of the time. 

The 3-ball tref oil_arc is obtained from the 3-sphere trefoil by deleting the star of a 
vertex. It contains the trefoil knot as a spanning arc and has optimal discrete Morse vector 
(1,0,0,0). The deletion of the star of a vertex from the 3-sphere double_tref oil yields the 
3-ball double_tref oil_arc with the double trefoil knot as spanning arc and optimal discrete 
Morse vector (1, 1, 1, 0). 

For the triple trefoil knot the deletion of a vertex from the 3-sphere triple_tref oil yields 
the 3-ball triple_tref oil_arc for which the optimal discrete Morse vector is (1,2,2,0); see 
Theorem El We found this vector in about 60% of the runs. 

A larger 3-ball knot that has the trefoil knot as spanning arc was constructed (via a pile 
of cubes) by Hachimori |32] , The best discrete Morse vector we found for knot is (1, 1, 1,0). 
It might as well be that knot admits (1,0,0,0) as optimal discrete Morse vector. The non- 
constructible 3-sphere nc_sphere [32j is obtained from knot by adding the cone over the 
boundary of knot. For this example, we found (1,0,0, 1) as optimal discrete Morse vector, 
but only in 12 out of 10000 runs. 

The triangulation bing is a 3-dimensional thickening of Bing's house with two rooms [13J 
due to Hachimori [32J (again, via a pile of cubes). It is a 3-ball with 480 vertices for which we 
found (1,0,0,0) as optimal discrete Morse vector in only 7 out of 10000 runs. We therefore 
can regard this ball as barely collapsible. 

A non-(4, 2)-colorable triangulation non_4_2_colorable of the 3-sphere was constructed 
in [2] with 167 vertices by using 10 copies of the double trefoil knot. The best discrete Morse 
vector we found once for this example in 10000 runs is (1,2,2,1). The average number of 
critical cells for non_4_2_colorable computed and normalized over only 10 random runs (for 
the sake of simplicity) as listed in Table [3] is roughly 25.2. 

5.6 Barycentric subdivisions 

Examples: trefoil_bsd, double_tref oil_bsd, triple_tref oil_bsd. 

Interestingly, the barycentric subdivisions trefoil_bsd, double_tref oil_bsd, and 
triple_tref oil_bsd of the knotted spheres trefoil, double_tref oil, and triple_tref oil, 
respectively, have a lower observed spectrum than the corresponding original spheres; com- 
pare Table [2j 
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5.7 Standard and exotic PL structures on 4-manifolds 

Examples: CP2, RP4, K3_16, K3_17, RP4_K3_17, RP4_llS2xS2. 

Freedman's classification [29] of simply connected closed topological 4-manifolds settled the 4- 
dimensional topological Poincare conjecture. The 4-dimensional smooth Poincare conjecture, 
however, is still wide open: Does the 4-dimensional sphere S have a unique differentiable 
structure or are there exotic 4-spheres that are homeomorphic but not diffeomorphic to 5 4 ? 
The categories PL and DIFF coincide in dimension 4 (see the survey of Milnor [61] and the 
references contained therein), and the 4-dimensional smooth Poincare conjecture therefore can 
be rephrased on the level of triangulations: Is every triangulation of S 4 PL homeomorphic 
to the boundary of the 5-simplex? 

Exotic structures on simply connected 4-manifolds have been intensively studied over 
the past years. One main task has been to find smaller and smaller k and I such that the 
connected sum (#/cCP 2 )#(-#ZCP 2 ) has exotic structures. While it is now known that 
CP 2 #(— #2 CP 2 ) [4J admits (infinitely many) exotic structures, the remaining interesting 
open cases are CP 2 #(— CP), CP 2 , and S 4 (the smooth Poincare conjecture). 

The example CP2 in Table [3] is the unique vertex-minimal 9- vertex triangulation of CP 2 
due to Kiihnel and Banchoff [42J and it carries the standard PL structure. 

The constructions of exotic structures are often delicate and it is not straightforward to 
derive corresponding triangulations. A very explicit example, though not simply connected, 
is due to Kreck [30] . 

Theorem 17 (Kreck |4QJ) The A- dimensional manifolds MP 4 #K3 and RP 4 #(S 2 x S 2 )* 11 
are homeomorphic but not diffeomorphic; the constituting components being equipped with the 
standard smooth structures. 

A 17-vertex triangulation K3_17 of the K3 surface with the standard PL type is due to 
Spreer and Kiihnel [67] . A vertex-minimal 16- vertex triangulation K3_16 of the topological 
K3 surface was previously found by Casella and Kiihnel [18]. It is not clear whether the two 
triangulations are PL homeomorphic — we tried bistellar flips to establish a PL homeomor- 
phism between these two triangulations, but without success. 

A vertex- minimal 16- vertex triangulation RP4 of RP 4 was obtained in [55] by applying 
bistellar flips to the 31-vertex standard triangulation of MP 4 by Kiihnel [41] . 

Let K and L be two triangulated 4-manifolds K and L with n and m vertices, respectively. 
Their connected sum K#L (or K# — L in cases when orientation of the components matters) 
is obtained from K and L by removing a 4-simplex from each of the triangulations and then 
gluing together the remainders along the respective boundaries. The resulting triangulation 
K#L then has n + m — 5 vertices. Triangulations of connected sums (S 2 x S 2 )* k , k>2, arc 
therefore easily constructed from a vertex-minimal 11- vertex triangulation of S 2 x S 2 [55] by 
taking connected sums and then applying bistellar flips to reduce the numbers of vertices. 
This way, we obtained triangulations of(S 2 xS 2 )^ 2 with 12 vertices (vertex-minimal; c.f. [55J), 
of (S 2 x S 2 )# 3 with 14 vertices, of (S 2 x S 2 )# 5 with 16 vertices, of (S 2 x S 2 )* 6 with 16 vertices, 
of (S 2 x S 2 )* 9 with 18 vertices, and of (S 2 x S 2 )* 11 with 20 vertices. 

Theorem 18 Let MP 4 , K3, and (S 2 x S 2 )* 11 be equipped with their standard PL structures. 
The PL A-manifold WP 4 #K3 has a triangulation RP^_K3_17 with 16+17—5 = 28 vertices and 
the PL 4-manifold MP 4 #(S 2 x S" 2 )* 11 has a triangulation RP4_llS2xS2 with 16 + 20-5 = 31 
vertices. While the underlying topological manifolds of these PL manifolds are homeomorphic, 
the respective triangulations are not PL homeomorphic. 
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By Theorem [18] we see that homeomorphic but not PL homeomorphic triangulations of 
4-manifolds can be constructed with only few vertices. (Most likely, the explicit numbers 
of vertices in Theorem [18] can be further reduced with bistellar flips. However, this would 
require a rather extensive search, which is beyond the scope of this article.) 

Theorem 19 The examples CP2, RP4, K3A6, K3_17, RP4_K3_17, and RP4_llS2xS2 have per- 
fect discrete Morse functions with 3, 5, 24, 24, 27, and 27 critical cells, respectively. 

Interestingly, the computed discrete Morse spectra of K3_16 and K3_17 look rather similar. 
The same can be said for the pair RP4_K3_17 and RP4_llS2xS2. 

5.8 Horn complexes 

Examples: Hom_C5_K4, Hom_n9_655_compl_K4, Hom_C6_compl_K5_small , Hom_C6_compl_K5 , 
Hom_C5_K5 . 

Horn complexes of certain graphs provide interesting examples of prodsimplicial manifolds |23j. 
The prodsimplicial structure allows to easily triangulate these manifolds without adding new 
vertices. 

The 3-dimensional Horn complex Hom_C5_K4 is a triangulation of the 3-dimensional real 
projective space RP 3 , the Horn complex Hom_n9_655_compl_K4 triangulates (S 2 x S 1 )* 13 . 

The 4-dimensional example Hom_C6_compl_K5_small with / = (33, 379, 1786, 2300, 920) is 
obtained from Hom_C6_compl_K5 with / = (1920, 30780, 104520, 126000, 50400) via bistellar 
flips. Both examples triangulate (5 2 x5 2 )* 29 , the first with computed normalized average 
63.92, the latter with normalized average 83.0. In only three out of 2000 runs we found the dis- 
crete Morse vector (1, 1, 59, 0, 1) for Hom_C6_compl_K5, but never the optimum (1, 0, 58, 0, 1). 
In contrast, both the lex and the rev_lex heuristics yielded (1,0,58,0, 1). In order to keep 
the list short, Table [3] only lists 10 random runs for Hom_C6_compl_K5. 

The Horn complex Hom_C5_K5 with / = (1020, 25770, 143900, 307950, 283200, 94400) is a 
triangulation of S" 3 x S 2 with normalized average 4.6. 

5.9 Higher-dimensional manifolds 

Examples: poincare, hyperbolic_dodecahedral_space , S2xpoincare, SU2_S03, RP5_24, 
non_PL, _HP2. 

The 16-vertex triangulation poincare |14] of the Poincare homology 3-sphere with /-vector 
/ = (16, 106, 180, 90) has the binary icosahedral group as its fundamental group. Since this 
group is non-cyclic, we have C2 > 2, and therefore every discrete Morse vector for poincare 
must have at least six critical cells, with (1,2,2, 1) being the optimal discrete Morse vector 
according to Table [3] cf. also Lewiner [45]. 

For the 21-vertex triangulation hyperbolic_dodecahedral_space [60 of the Weber- 
Seifert hyperbolic dodecahedral space [71] with face vector / = (21,193,344,172) the best 
discrete Morse vector we found is (1, 4, 4, 1). The fundamental group of this manifold can be 
presented with 4 generators; see Table [5] 

The product triangulation S2xpoincare of S 2 (taken as the boundary of a tetrahedron) 
with poincare again has the binary icosahedral group as its fundamental group, inherited 
from poincare; for constructing product triangulations see [54] and references therein. The 
best discrete Morse vector we found for this examples in 103 out of 1000 runs is (1, 2, 3, 3, 2, 1). 
(Table [3] list only 20 random runs for S2xpoincare to keep the list short.) 
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The two 5-manifolds SU2/S03 and MP 5 have homology vectors (Z, 0, Z, Z, 0, Z) and 
(Z, Z2, 0, Z2, 0, Z) and triangulations SU2_S03 with 13 vertices and RP5_24 with 24 vertices, 
respectively [49\ I55j. The 15- vertex triangulation _HP2 of an 8-dimensional manifold 'like a 
quaternionic projective plane' by Brehm and Kiihnel [TE] has homology (Z, 0, 0, 0, Z, 0, 0, 0, Z). 

Theorem 20 The triangulations SU2J503, RP5_2^., and _HP2 have optimal discrete Morse 
vectors (1, 0, 1, 1, 0, 1), (1,1,1,1,1,1), and (1,0,0,0,1,0,0,0,1), respectively. 

The 18-vertex non-PL triangulation non_PL [14J of the 5-dimensional sphere S 5 admits 
(1, 0, 0, 2, 2, 1) as discrete Morse vector. 



5.10 Random 2-complexes and fundamental groups 
Example : rand2_n25_p0 . 328 

In generalization of the classical Erdos-Renyi model for random graphs [26], Linial and 



Meshulam [37] considered random 2-dimensional complexes with complete 1-skeleton on 
n vertices; every triangle with vertices from the set {1, ... , n} is then added with proba- 
bility p independently. Let Y(n,p) be the set of such complexes. For the elements of Y(n,p), 
Linial and Meshulam proved a sharp threshold for the vanishing of the first homology with 
Z2-coefficients, 



lim Prob[ Y e Y (n, p) \ H x (Y, Z 2 ) = ] = < 



1 for p = 21ogn+a,(n) 



for p 



21ogra— ui(n) 



for any function u(n) — > 00 as n — > 00. Replacing homological connectivity by simple 
connectivity, Babson, Hoffman, and Kahle [6] showed that there is a range for p for which 
asymptotically almost surely the complexes Y £ Y(n,p) have non-trivial fundamental groups 
with trivial abelianizations, 

lim ProbfY €Y(n,p)\ir 1 (Y) = Q] =1 for p> ( gggjgfch s 

with the exponent \ being best possible. 

More recently, Cohen et al. [21] showed that for p <C n^ 1 asymptotically almost surely 
the complexes Y G Y(n,p) admit a discrete Morse functions with no critical 2-cells. 

The example rand2_n25_p0 . 328 on n = 25 vertices from Table[3]with homology (Z, 0, Z 475 ) 
has 751 triangles, each picked with probability p = 0.328. We found the optimal discrete 
Morse vector (1,0,475) in 275 out of 10000 runs. 

According to Seifert and Threlfall [661 §44], a presentation of the fundamental group of 
a simplicial complex can be obtained via the edge-path group. For this, a spanning tree 
of edges is deleted from the 1-skeleton of the complex and each remaining edge contributes 
a generator to the fundamental group, while each triangle of the 2-skeleton contributes a 
relator; see [58] for an implementation. 

We used the GAP command Simplif iedFpGroup to simplify the edge-path group presen- 
tation of the fundamental group. The heuristic Simplif iedFpGroup not necessarily outputs a 
minimal presentation of a finitely presented group with the minimal number of generators and 
relators. Nevertheless, even in the case of huge complexes, Simplif iedFpGroup succeeded 
with recognizing trivial, cyclic (one generator, at most one relator), or free groups (no rela- 
tors). In Table [U we list for the examples of Table [3] the number of generators (Ge.) and the 
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number of relators (Re.) of the initial presentation and the number of generators (SGe.) and 
the number of relators (SRe.) of the simplified group along with the resulting fundamental 
group (F. Gr.) and the time it took for the simplification. In Tables [5H3 F(k) denotes the 
free group with k generators. 

In the Tables [6] and [7] we list resulting fundamental groups for random 2-complexes with 25 
and 50 vertices, respectively. In these tables, the Linial-Meshulam threshold can be observed 
quite nicely. For p = ^f 25 ~ 0.26 and p = 21 g| 50 « 0.16, 73 and 75 out of 100 random 
examples with n = 25 and n = 50 vertices had trivial fundamental groups, respectively. 
Thus, for these values of p we precisely are in the range of the slope of the threshold. While 
most of the examples in the Tables and [7] have free fundamental groups, we found 'non- 
free' examples (for which their presentations could not be simplified to remove all relators) 
in the range when p is slightly smaller than -, the value for which Linial, Meshulam, and 
Rosenthal [H] constructed acyclic examples as sum complexes. 

In our experiments we did not observe the Babson-Hoffmann-Kahle examples with non- 
trivial fundamental groups that have trivial abelianizations. However, as pointed out by 
Kenyon, 'exceptional events' can occur for random groups in the case when n is small while 
the asymptotical behavior can be rather different; cf. [62} pp. 42-43]. 

5.11 Vertex-homogeneous complexes and the Evasiveness Conjecture 
Example : contractible_vertex_homogeneous . 

As remarked by Kahn, Saks, and Sturtevant [36] we have the following implications for 
simplicial complexes: 

non-evasive ==>• collapsible ==> contractible ==> Z-acyclic. 

The Evasiveness Conjecture [36J for simplicial complexes states that every vertex-homoge- 
neous non-evasive simplicial complex is a simplex. The first examples of vertex-homogeneous 
Z-acyclic simplicial complexes different from simplices were given by Oliver (cf. [36]); see [50] 
for further examples. While join products and other constructions can be used to derive 
vertex-homogeneous contractible simplicial complexes different from simplices, non-trivial 
vertex-homogeneous non-evasive examples cannot be obtained this way [72J. 

The smallest example contract ible_vertex_homogeneous of a contractible vertex-homo- 
geneous simplicial complex from |50] is 11-dimensional with 

/ = (60, 1290, 12380, 58935, 148092, 220840, 211740, 136155, 59160, 16866, 2880, 225). 

The best discrete Morse vector we found with the lex and the rev_lex heuristics for this 
contractible space is (1,0,0,4,8,4,0,0,0,0,0,0). We do not know whether the example is 
collapsible or not. 

Acknowledgment s 

Thanks to Karim Adiprasito, Herbert Edelsbrunner, Alex Engstrbm, Michael Joswig, Roy 
Meshulam, Konstantin Mischaikow and Vidit Nanda for helpful discussions and remarks. 



17 



Figure 4: The graph A^+q with k + 6 vertices. 

Appendix A: Complexes on which our heuristic fails 

In this section, we construct simplicial complexes on which our random approach will most 
likely go far from guessing the right Morse vector. On these examples, exponentially many 
rounds (in the number of facets) of the program may be necessary, before an optimal Morse 
vector shows up as output. Such pathological examples can be produced in any positive 
dimension. The crucial idea is highlighted by the following 1-dimensional case. 

Example 21 Let k be a positive integer. Let A^+q be the graph consisting of two cycles 
of length 3 that are connected by a path of k edges; see Figure 01 Since A^+e has no free 
edge, our algorithm picks an edge e uniformly at random and removes it. The final outcome 
depends on this choice, and on this choice only: 

• If e belongs to the A:-edge path, it is easy to see that the program will always output 
the discrete Morse vector (2,3). 

• If instead e belongs to one of the two triangles, then the program will always output 
the Morse vector (1,2). 

Hence, the algorithm finds a perfect Morse function on ^4^+6 with probability p = qt^- For 
large k, the algorithm will most likely (i.e. with probability q = gr^) return a Morse vector 
that is 'off by 2', displaying 5 critical cells instead of 3. 

Example 22 Let s be a positive integer. Let Bk +e (s) be a bouquet of s copies of A^q. 
An optimal discrete Morse function on Bk+e(s) has Morse vector (1,2s). Finding a discrete 
Morse function on Bk + Q(s) is the same as (independently) finding s discrete Morse functions 
on the s copies of Ak+e- Therefore, the probability of getting the optimal Morse vector 
on 5fc + 6(s) is p s , where p = qt^- This corresponds to putting together s optimal Morse 
functions on the different copies of A^+e, or in other words, to picking one favorable edge in 
each copy of ^4^+6 • For < i < s, the probability that the program outputs the Morse vector 
(1 + i, 2s + i) is (^jp s ~ l (l — p) 1 , corresponding to i 'bad choices' and s — i 'good choices'. 

To show that an analogous phenomenon occurs also in higher dimensions, let us recall a 
classical definition in PL topology. 

Definition 23 Let C be a d-dimensional complex. A stacking operation on C is the transition 
from C to C = (C — star (a, C) ) U a * link(<7, C), where a is an arbitrary facet of C and a 
is a new vertex (e.g., the barycenter of a). More generally, we say that C is obtained from 
C by stacking if some finite sequence of stacking operations leads from C to C . 

Each stacking operation adds d facets; so, a complex obtained by performing s stacking 
operations on a d-simplex has exactly ds + 1 facets. 

Lemma 24 If C is obtained from a simplex by stacking, then C is shellable. In particular, 
it is endo- collapsible: For any facet a of C , there is a sequence of elementary collapses that 
reduces C — a to dC . 
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In dimension d > 3, there is no guarantee that any sequence of elementary collapses 
on C — a can be continued until one reaches dC . This is why, in the following example, 
probabilities have to be estimated rather than computed. 

Example 25 Let d, k be positive integers, with k = 1 mod d. Take a disjoint union of d + 1 
edges eo, • • • , &di and a (f-simplex S 1 (with facets Fq, . . . , F^). For each i in {0, ... , d}, glue 
in the boundary of the join i*^ * e\. The resulting complex is homotopy equivalent to a 
bouquet of d + 1 spheres of dimension d; the homotopy is just the contraction of the central 
simplex S to a point. Let be a complex obtained by stacking the simplex S exactly 
s times, so that S gets subdivided into k = ds + 1 simplices of the same dimension. 

Note first of all that coincides with the A^+q of Example ETJ Since Cf. has no free 
(d — l)-faces, our algorithm starts by removing some d-face a at random. We have two 
possible cases: 

• With probability ( d+2 )(d+i)+k we P^ a f rom the subdivision of the central simplex. 

• With probability ( rf+^Xd+Tj+fc we P^ ° f rom one °f the ^-spheres. 

In the first case, some sequence of elementary collapses reduces Cf — a onto d+ 1 disjoint 
d-spheres. So our algorithm will output a Morse vector that is either (1, 0, . . . , 1, d + 2), or a 
larger vector; but certainly not the vector (1, 0, . . . , 0, d + 1). 

Thus the probability of obtaining the optimal Morse vector (1, 0, . . . , 0, d + 1) is positive, 

but smaller or equal to ^t^{dt^)Xk • ^ s ^ & e ^ s l ar g er > this upper bound gets smaller and 
smaller. 

Example 26 By taking a bouquet of w copies of Example 1251 we obtain a complex Bf(w). 
For d = 1, Bf{w) coincides with the Bk+e(w) of Example 1221 The probability of seeing the 

perfect Morse vector (1,0,..., 0, (d + l)w) on B^(w) is smaller or equal to ( (d+2)(d+i)+k ) • 

For practical purposes, it is useful to understand how this probability grows with respect 
to the number N of facets. In fact, given a complex with TV facets, we would like to concretely 
know how often we should run the algorithm before we can expect an optimal Morse vector 
to appear among the outputs. 

For the sake of brevity, we do the calculations in dimension one — but similar estimates 
can be easily derived in all dimensions. The graph constructed in Example [22] has N = 
(6 + k)w edges. To study the probability (g^) w of finding an optimal Morse function, we 
should regard N as a constant, write w as gx£, an d study the function 



P(k) 



N 

6 \ 6+k 



6 + k 



Now, classical calculus reveals that the function x \ — > x x = e xlogx is strictly decreasing on 
the interval (0,e -1 ) and strictly increasing on (e _1 ,oo). It achieves its minimum at e . So, 
given any bijection g : (0, oo) — > (0,1), the function y i — > g(y) 9 ^ v ' achieves its minimum at 
the (unique) point y such that g{y) = e -1 . Applying this to g(y) = we get 

N JV 

Yet we wanted to minimize the function P(k) over the integers, not over the reals. However, 
if we choose the integer k so that is close to e _1 , one can see that the value of P(k) is 
close to P(e~ 1 ). The minimum is in fact achieved at k = 10. Thus P(k) can be as small as 
e~ cN , where c is some constant 'close' to g-r: It is in fact c = 4; (log 8 — log 3) sa 0, 0613018. 



mm 
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Appendix B: Library and Tables 



Table lists computational results for the examples of Section 

Of each example we present the discrete Morse spectrum we experimentally observed 
in a certain number of runs (usually 10000, when not otherwise stated; sometimes we did 
fewer runs for reasons concerning either excessive computation time or excessive variance of 
the spectrum). To sum up in a single coefficient the complicatedness of the spectrum, we 
introduce the following averages: 

• , the average number of critical cells for the vectors in the normalized discrete Morse 
spectrum of a simplicial complex C; 

• c~, the average number of critical cells for a collection of normalized discrete Morse 
vectors for C. 

These two coefficients are of some interest if we want to reduce the size of the complex as 
a preprocessing step for homology computations. In fact, the second value gives an estimate 
for the number of cells that we are left with for the Smith normal form computations. 

The longer we run Random Discrete Morse, the better the approximation of by 
c~ will get. If we had an infinite amount of time (and patience), the optimal discrete Morse 
vectors would eventually show up — as they are always contained in the respective (finite!) 
discrete Morse spectra. 

In Table [3j optimal discrete Morse vectors are highlighted in bold. We wrote an output 
vector in italics if it is the best we could find with our algorithm and we do not know if it is 
indeed the optimal or not. 

For Table [H we replaced the random choices in our algorithm with a deterministic lex- 
icographic or reverse lexicographic choice. The labeling of the vertices of course now plays 
a role; see |59] for a discussion of a randomized version (by randomly renumbering vertices 
first) of rev_lex. 

All computations were run on a cluster of 2.6 GHz processors. 



Table 3: Library of triangulations and discrete Morse spectra. 



Name of example/Homology/ 
/-vector/dY 



Distribution of obtained 
discrete Morse vectors 
in 10000 rounds 



Time for Hasse diagram/ 

Time per round 

(in Hour:Min:Sec.Frac) 




dunceJiat 



(1,1,1): 



10000 



0.004 
0.00024 




14.0000 

regular _2_21_23_1 




9722 
277 
1 








004 
00076 





9337 
649 
14 








008 
00201 





2185 
1874 
1847 
1265 
1048 
704 
318 
275 
140 








228 
00428 
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Table 3: Library of triangulations and discrete Morse spectra (continued). 



Name of examplc/Homology/ 
/- vector /c~ 



Distribution of obtained 
discrete Morse vectors 
in 10000 rounds 



Time for Hasse diagram/ 

Time per round 

(in Hour:Min:Sec.Frac) 



dunce_tiat_in_3_ball 

(Z, 0,0.0) 

(8,25,30,12) 

1.0000 

Barnette_sphere (non-polytopal) 

(Z,0,0.Z) 

(8,27,38,19) 

2.0000 

B_3_9_18 (non-shellable ball) 
(Z, 0,0,0) 
9.33,43,18 
1.0000 

tref oil_arc 

(Z, 0,0,0) 

(12,58,85,38) 

1.0952 

trefoil 

(Z,0,0,Z) 

(13,69,112,56) 

2.0778 

rudin (Rudin's ball) 
(Z, 0,0,0) 
(14,66,94,41) 
1.0000 

double_tref oil_arc 

(Z, 0,0,0) 
^3,145,66) 



Domcare 



(Z,0,0.Z) 

(16,106,180,90) 

6.1952 



double_tref oil 

Z,0,0,Z) 

16,108,184,92) 

.5338 



2,4, 
2,5, 
(2,6, 



478) 
479 
480 



(2,7,481' 
1,9,484) 



477 
8,482 
10,48 
9,483): 
11,486* 
10,484* 
11,485* 
6,479):' 
8,481): 



(1,0,0,0): 
(1,0,0,1): 
(1,0,0,0): 



1,0,0,0) 

1,1,1,0 

1,2,2,0); 

1,0,0,1; 

1,1,1,1 

l,2,2,lj 
(1,0,0,0): 



1,1,1,0, 

(1,2,2,0 



1,3,3 
2,3,2 
1,4,4 
(2,4,3 

(1,2,2 

(1,3,3 
1,4,4 
2,4,3 
(2,3,2 
(1,5,5 

(1,1,1 
1,0,0 

1,2,2 
1,3,3 

(1,4,4 
2,3,2 

(2,4,3 



66 
66 
54 
41 
40 
24 
21 
12 
8 
4 
3 
3 
1 
1 

10000 



9529 
466 
5 

9617 
377 
6 

10000 



7080 
2698 
197 
18 
6 
1 

9073 
864 
45 
7 
6 
5 

4550 
3972 
1316 
145 
8 

7 

9 



0.004 
0.00049 



10000 0.004 

0.00060 



10000 0.004 

0.00073 



0.004 
0.00158 



0.004 
0.00208 



0.004 
0.00107 



0.012 
0.00329 



0.016 
0.00400 



0.012 
0.00408 
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Table 3: Library of triangulations and discrete Morse spectra (continued). 



Name of examplc/Homology/ 
/- vector /c~ 



Distribution of obtained 
discrete Morse vectors 
in 10000 rounds 



Time for Hasse diagram/ 

Time per round 

(in Hour:Min:Sec.Frac) 



triple_tref oil_arc 

(Z, 0,0,0) 

(17,127,208,97) 

5.9352 



triple_tref oil 

Z,0,0,Z) 

18. 143,250,125) 

.9^98 



hyperbolic_dodecahedral_space 

(2,Z3,0,Z) 
21,190,338,169 
11.4672 



S_3_50_1033 (random) 
(Z,0,0,Z) 

50,1083,2066,1033) 
3.1966 



S_3_100_4850 (cyclic polytope) 
(Z,0,0,Z) 

100,4950,9700,4850) 
2.0000 

600_cell 

(Z,0,0,Z) 
120,720,1200,600) 
2.0000 

non_4_2_colorable 

(Z,0,0,Z) 

(167,1579,2824,1412) 
25.2 



(1,2,2 



3,3 
4,4 



5,5 

4 3 
3 2 

5 4 
6,6 

2,2 

1,1 

3,3 



.: 

6 6 



its 

6,6 
7,7 
5 4 
6,5 

8 8 
7,6 
8,7 

9 9 
9,8,., 
10, lo. 



1): 



10 1) 



(1,0,0,1): 



(1,0,0,1): 



4,15,12 1): 



(7,20 



1 



11,1 
12,1 
11 

13; 1 
10 1 

14: 1 

14 



6027 
3220 
569 
77 
51 
42 
10 
4 

4427 
3080 
1911 
430 
57 
40 
33 
15 
4 
3 

4792 
3338 
1245 
326 
82 
80 
62 
45 
18 
8 
3 
1 

7087 
1383 
697 
386 
189 
118 
42 
25 
18 
14 
12 
9 
9 
4 
3 
2 
1 
1 

10000 



10000 



0.024 
0.00528 



0.024 
0.00640 



0.036 
0.01017 



0.900 
0.153 



17.829 
1.883 



0.364 
0.076 



1.728 
0.254 

[10 rounds] 
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Table 3: Library of triangulations and discrete Morse spectra (continued). 



Name of examplc/Homology/ 
/- vector /c~ 



Distribution of obtained 
discrete Morse vectors 
in 10000 rounds 



Time for Hasse diagram/ 

Time per round 

(in Hour:Min:Sec.Frac) 



Hom_C5_K4 (RP 3 ) 

(Z,Z 2 ,0,Z) 

(240, 1680, 2880, 1440) 

4.0496 

tref oil_bsd 

(Z,0,0,Z) 

(250,1594,2688,1344) 
2.0202 

knot 

(Z, 0,0,0) 

(380.1929,2722,1172) 
3.1194 



nc_sphere 

(Z,0,0,Z) 

^381^2309, 3856, 1928) 



double_tref oil_bsd 

(Z,0,0,Z) 

^400,2608,4416,2208) 



£ in j? 

(Z, 0,0,0) 

(480,2511,3586,1554) 
3.0456 

triple_tref oil_bsd 

(Z,0,0,Z) 

(536,3536,6000,3000) 
5.7352 



S_3_1000_2990 (stacked sphere) 
(Z,0,0,Z) 

1000,3990,5980,2990) 
2.0000 

Hom_n9_655_compl_K4 ((S 2 x S 1 )* 13 ) 
(Z,Z 13 ,Z 13 ,Z) ' ' 

(3096,22104,38016,19008) 
28.68 



CP2 

(Z,0,Z,0,Z) 

(9,36,84,90,36) 

3.0012 



1,1,1,1. 

1,2,2,1 
3,2,1 
3,3,1 



\\ 

1 
1 
1 

2 
2 
1 
2 
(1 

(1 



0,0,1, 

1,1,1, 
2,2,1 



1,1, Oh 
2,2,0 : 



3,2 
3,3 
4,3 

1,1 
2,2 
3,3 
4,4 
0,0 
3,2 
6,6 
4,3 
5,4 
5,5 

1,1 
0,0 

2,2 
3,3 
4,4 
3,2 
4,3 

1,1 
2,2 
0,0 

3,3 
3,2 

2,2 

1,1, 
3,3 
4,4 
5,5 
3,2 
4,3 
6,6 
5,4 
8,8 

0,0 



: 

: 

0) : 

1) : 

1 : 
1 : 
1): 

1): 

1): 
1 : 
1 : 
1 : 
1): 

1): 
1): 
1 : 
1 : 
1 : 
1 : 
1): 

: 

: 
0): 

0): 

0) : 

% 

1 : 
1 : 
1 : 
1 : 
1 : 
1 : 

1) : 

1): 



(1,13,13,1) 

(1,14,14, r 

(1,15,15,1 
(2,14,13,1 
(2,15,14,1 
(2,16,15,1 

(1,0,1,0,1) 

(1,1,2,0,1) 



9753 
240 
6 
1 

9902 
95 
3 

9414 
560 
15 
9 
2 

7902 
1809 
234 
25 
12 
9 
3 
3 
2 
1 

4819 
4274 
833 
64 
4 
4 
2 

9764 
217 
7 
6 
6 

4793 
3390 
1543 
208 
22 
20 
17 
3 
3 
1 

10000 



67 
20 
5 
5 
2 
1 

9994 
6 



1.864 
0.379 



1.716 
0.308 



1.576 
0.813 



3.228 
0.470 



4.376 
0.811 



2.788 
1.398 



8.024 
1.456 



8.444 
1.498 



5:39.809 
45.682 
[100 rounds] 



0.012 
0.00226 
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Table 3: Library of triangulations and discrete Morse spectra (continued). 



Name of examplc/Homology/ 
/- vector /c~ 



Distribution of obtained 
discrete Morse vectors 
in 10000 rounds 



Time for Hasse diagram/ 

Time per round 

(in Hour:Min:Sec.Frac) 



RP4 

(Z,Z 2 ,0,Z 2 ,0) 
16,120,330,375,150) 
5.0490 



K3_16 (unknown PL type) 
(Z,0,Z 22 ,0,Z) 
16,120,560,720,288) 
24.8218 



K3_17 (standard PL type) 
(Z,0,Z 22 ,0,Z) 
17,135,610,780,312) 
24.8978 



RP4_K3_17 

(Z,Z 2 ,Z 22 ,Z 2 ,0) 

(28,245,930,1150,460) 

28.56 



RP4_llS2xS2 

(Z,Z 2 ,Z 22 ,Z 2 ,0) 
31,283,1052,1295,518) 
28.46 



1,1 



1,1,1. 

2,1,1. 
2,2,1, 
3,1,1 
3,2,1, 
3,3,1 
2,1,1 

22,0 

23,0 
24,0 
25,0 
23,1 
24,1 
24,2 
25,3 
25,1 
24,0 
26,4 
26,3 
26,0 
27,5 
27 
27 
26 
28 
28 
27 
29 
26 



(1 
1 
1 

1 
1 

(2 

(1 

(1 
1 
1 
1 
1 
1 

1 
1 

2 

1 
1 

1 
1 
1 
1 

1 
1 

1 
1 
1 

(1 
(1 

fi 
i 
i 

i 
i 

(1,0 

2,3 
fl,0 

1,0 

1,1 

1,2 
(1,5 

(1,1 

1,2 
1,3 

1,1 
1,2 
1,4 

1, 1 
1,3 

1,3 
1,3 
1,5 
(2,4 

(1,1 

1,2 
1,3. 

1,1,24,2 
(1,1,25,3 
(1,1,26,4 
1,2,27,4 
1,4,26,1 
2,4,25,1 



22,0 

23,0 
24,0 
25,0 
26,0 
23,1 
25,3 
24,0 
24,2 
26,4 
24,1 
27,3 
27,0 

23,1 

24,1 
25,1 
24,2 
25,2 
26,1 
26,4 
26,2 
27,3 
28,4 
30,4 
25,1 

23,1 

24,1 
25,1 



9765 
136 
89 
5 
3 
1 
1 

6702 
2615 
506 
60 
31 
15 
13 
9 
6 
5 
4 
4 
4 
3 
3 
2 
2 
2 
1 
1 
1 
1 

6337 
2939 
618 
78 
8 
6 
4 
3 
2 
2 
1 
1 
1 

55 
24 
8 
3 
2 
2 
1 
1 
1 
1 
1 
1 

51 
29 
14 
1 
1 
1 
1 
1 
1 



0.056 
0.01678 



0.168 
0.04417 



0.196 
0.05093 



0.392 
0.0754 
[100 rounds] 



0.496 
0.0945 
[100 rounds] 
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Table 3: Library of triangulations and discrete Morse spectra (continued). 



Name of example/Homology/ 
/-vector/dY 



Distribution of obtained 
discrete Morse vectors 
in 10000 rounds 



Time for Hasse diagram/ 

Time per round 

(in Hour :Min: Sec. Frac) 



Hom_C6_compl_K5_small ((S 2 x S 2 )* 29 ) 

33,379,1786,2300,920) 
'3.92 



Hom_C6_compl_K5 ((S 2 X S 2 )* 29 ) 
(Z,0,Z 58 ,0,Z) 

(1920, 30780, 104520, 126000, 50400) 
83.0 



SU2_S03 

(Z, 0, Z, Z, 0, Z) 

(13,78,286,533,468,156) 

4.1354 



non_PL 

(Z,0,0,0,0,Z) 

(18.139,503,904,783,261) 

6.1328 



RP5_24 

Z,Z 2 ,0,Z 2 ,0,Z) 
24,273,1174,2277,2028,676) 
1766 



S2xpoincare 

(Z,0, Z, Z.0.Z) 

64, 1156, 5784, 11892, 10800, 3600) 
15.70 



S_5_100_472 (stacked sphere) 
(Z,0,0,0,0,Z) 

(100, 579, 1430, 1895, 1416, 472) 
2.0000 



1,59,0,1 
2,60,0,b 
3,61,0,1' 
0,58,0,1 



62,0, 
63,0, 



6,64,0. 



60,0. 

63, i: 



7,65,0 : 



10,68,0,1 
17,75,0,1 
7,65,0,1" 
8,66,0,1 
9,67,0,1 
11,69,0,1 
16,73,0,1 



(1,0 



0,2,2.. 
' 3 



3,3,2 : 
' 4 



1 

1' 

1' 
1' 
1' 
1' 

r 
r 
r 

i 

r 
r 
i' 
r 
r 
r 
r 
r 
r 
r 
r 
r 
i' 

i 
i 

"i 



33 
30 
12 
11 
5 
3 
3 
1 
1 
1 

3 
2 
1 
1 
1 
1 
1 

9369 
554 
35 
32 
5 
4 
1 

9383 
441 
134 
23 
12 
2 
2 
1 
1 
1 

9181 
344 
315 
97 
21 
15 
9 
6 
5 
3 
1 
1 
1 
1 

6 
3 
2 
2 
2 
2 
1 
1 
1 

10000 



1.460 
0.348 

[100 rounds] 



2:18:33.603 
19:26.475 
[10 roundsl 
([2000 rounds] 



0.124 
0.03250 



0.324 
0.06964 



Sec.l5~B 



1.800 
0.429 



46.611 
8.460 

[20 rounds] 
([1000 rounds], 



1.188 
0.309 



Sec. 153 
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Table 3: Library of triangulations and discrete Morse spectra (continued). 



Name of examplc/Homology/ 
/- vector /c~ 



Distribution of obtained 
discrete Morse vectors 
in 10000 rounds 



Time for Hasse diagram/ 

Time per round 

(in Hour:Min:Sec.Frac) 



Hom_C5_K5 (S 3 xS 2 ) 
(Z,0,Z, Z,0,Z) 

1020, 25770, 143900, 307950, 283200, 
94400) 

4.6 
_HP2 

(Z,0,0,0,Z,0,0,0,Z) 

(15, 105, 455, 1365, 3003, 4515, 4230, 

2205, 490) 
3.1212 



contract ible_vertex_tiomogeneous 

(Z, 0,0, 0,0, 0,0, 0,0, 0,0,0) 
(60, 1290, 12380, 58935, 148092, 

220840, 211740, 136155, 59160, 

16866,2880,225) 
273.6 



(1,0,1,1,0,1 

(1,0,2,2,0,1 
1,1,2,1,0,1 



(1,0,0,0, 

(1,0,0,1, 



0,0,2, 
0,0,3, 
0,0,0, 
0,1,2, 
0,1,4, 
0,1,1, 
0,0,1, 

0, 0, 1, 
0,0,4, 



1,0,0,0,1) 

2,0,0,0,1 
3,0,0,0,1 
4,0,0,0,1 
2,1,0,0,1 
2,0,0,0,1 
4,0,0,0,1 
1,0,0,0,1 
5,3,0,0,1 
1,0,0,0,1 
6,4,0,0,1 
5,0,0,0,1 



9474 
459 
46 
7 
4 
3 
2 
1 
1 
1 
1 
1 



3,38,98,83,20,0,... 

3,42,72,43,10,0,... 

4,34,80,64,14,0,... 

4,53,87,56,18,0,... 

5,63,110,61,9,0,... 

5,70,139,92,18,0,.. 

6,42,115,108,29,0,. 

8,75,160,113,20,0, 

9,66,124,89,22,6,.. 

13,74,144,97,14,0, 



,0 
,0 

,0 

:8 

.01 



: 1 
..0 : 1 

,6): 1 
.,6): 1 



16:16:14.156 
2:53:37.911 
[10 rounds] 



11.348 
1.425 



11:42:37.994 
1:39:02.745 
[10 rounds] 



Table 4: Discrete Morse vectors with lex and rev_lex heuristics. 



Name of example 



lex 



rcvJcx 



dunceJiat 
d2nl2g6 

regular _2_21_23_1 

rand2_ri25_p0 . 328 

dunce_tiat_in_3_ball 

Barnette_sphere 

B_3_9_18 (non-shellable ball) 

tref oil_arc 

trefoil 

rudin (Rudin's ball) 
double_tref oil_arc 
poincare 
double_tref oil 
triple_tref oil_arc 
triple_tref oil 

hyperbolic_dodecahedral_space 
S_3_50_1033 (random) 
S_3_100_4850 (cyclic polytope) 
600_cell 

non_4_2_colorable 
Hom_C5_K4 (RP 3 ) 
tref oil_bsd 
knot 

nc_sphere 

double_tref oil_bsd 
bing 

triple_tref oil_bsd 
S_3_1000_2990 (stacked sphere) 
Hom_n9_655_compl_K4 ((S 2 x S 1 )* 13 ) 



>) 



1,12,1 
(1,30,1) 

(1,0,476 
(1,0,0,0) 
1,0,0,1 
(1,0,0,0) 

(1,2,2,0) 
(1,2,2,1) 
(1,0,0,0) 
(1,3,3,0) 

1,2,2,1) 

1,3,3,1 

\mi 

1,0,0,1 
(1,0,0.1) 

(1,30,30,1) 

1,1,1,1) 

(1,2,2,1) 

1,2,2,1 
1,3,3,1 

1,1,1,0 
1,4,4,1 
(1,0.0,1) 

2,14,13,1) 



(1, 
(1, 
(1, 
(1, 
(1, 
(1, 
(1, 
(1, 
(1, 
fl, 

(1, 
(1, 
(1, 
1, 
1- 
(1, 
(1, 
(1, 
(1, 
fl, 

(1, 
(1, 

(1. 

1, 
1, 
1, 
(1- 

u, 

(1, 
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Table 4: Discrete Morse vectors with lex and rev_lex heuristics (continued). 



(1,0,1,0,1) (1,0,1,0,1) 

1,1,1,1,1), 1,1,1,1,1) 



Name of example lex revJex 

CP2 
RP4 

K3_16 (unknown PL type) (i;o;22,0,i) (1^0^22,0,1) 

K3_17 (standard PL type) (l,0,22,0,l) (1,0,22,0,1) 

RP4_K3_17 (1,1,23,1,1) (1, 1,23,1,1 

RP4_llS2xS2 (1,1,23,1,1) (1,1,23,1,1) 

Hom_C6_compl_K5_small ((5 2 x5 2 ) #29 ) (l,0,58,0,l) (l,0,58,0,l 

Hom_C6_compl_K5 ((S 2 x S 2 )* 29 ) (l,0,58,0,l) (l,0,58,0,l) 

SU2_S03 (1,0,1,1,0,1) (1,0,1,1,0,1) 

non_PL (1,0,0,2,2,1) (1,0,0,2,2,1) 

RP5_24 (l, 1,1,1,1,1) (1,1,1,1,1,1) 

S2xpoincare (1,2,3,3,2,1) (1,2,3,3,2,1) 

S_5_100_472 (stacked sphere) (1,0,0,0,0,1) (1,0,0,0,0,1) 

Hom_C5_K5 (S* x S 2 ) (1,0,1,1,0,1) (1,0,1,1,0,1) 

-HP2 (1,0,0,0,1,0,0,0,1) (1,0,0,0,1,0,0,0,1) 

contractible_vertex_tiomogeneous (1,0,0,4,8,4,0,0,0,0,0,0) (1,0,0,4,8,4,0,0,0,0,0,0) 



Table 5: Simplified presentations of fundamental groups with GAP. 



Name of example 


Ge. 


Re. 


SGc. 


SRc. 


F. Gr. 


Time 


dunce_tiat 


17 


17 











0.036 


d2nl2g6 


55 


44 


12 


1 


12 gen. 


0.132 


regular _2_21_23_1 


127 


98 


30 


1 


30 gen. 


0.304 


rand2_ri25_p0 . 328 


276 


751 











0.876 


dunce_tiat_in_3_ball 


18 


30 











0.044 


Barnette_sphere 


20 


38 











0.048 


B_3_9_18 (non-shellable ball) 


25 


43 











0.048 


tref oil_arc 


47 


85 











0.092 


trefoil 


57 


112 











0.084 


rudin (Rudin's ball) 


53 


94 











0.060 


double_tref oil_arc 


79 


145 











0.148 


poincare 


91 


180 


2 


2 


2 gen. 


0.104 


double_tref oil 


93 


184 











0.160 


triple_tref oil_arc 


111 


208 











0.164 


triple_tref oil 


126 


250 











0.156 


hyperbolic_dodecariedral_space 


170 


338 


4 


5 


4 gen. 


0.276 


S_3_50_1033 (random) 


1034 


2066 











6.372 


S_3_100_4850 (cyclic polytope) 


4851 


9700 











2:38.918 


600_cell 


601 


1200 











2.220 


non_4_2_colorable 


1413 


2824 











12.644 


Hom_C5_K4 (RP 3 ) 


1441 


2880 


1 


1 


z 2 


13.492 


tref oil_bsd 


1345 


2688 











11.292 


knot 


1550 


2722 











12.536 


nc_spriere 


1929 


3856 











23.509 


double_tref oil_bsd 


2209 


4416 











30.002 


bing 


2032 


3586 











22.877 


triple_tref oil_bsd 


3001 


6000 











57.947 


S_3_1000_2990 (stacked sphere) 


2991 


5980 











58.099 


Hom_n9_655_compl_K4 ((S 2 X S 1 )* 13 ) 
CP2 


19009 


38016 


13 







59:14.006 


28 


84 








0.036 


RP4 


105 


330 


1 


1 


z 2 


0.200 


K3_16 (unknown PL type) 


105 


560 











0.344 


K3_17 (standard PL type) 


119 


610 











0.372 


RP4_K3_17 


218 


930 








z 2 


0.408 


RP4_llS2xS2 


253 


1052 








z 2 


0.472 


Hom_C6_compl_K5_small ((S 2 xS 2 )* 29 ) 


347 


1786 











2.420 


Hom_C6_complJ(5 ((S 2 xS 2 )* 29 ) 


28861 


104520 











3:57:09.873 


SU2_S03 


66 


286 











0.180 


non_PL 


122 


503 











0.368 


RP5_24 


250 


1174 


1 


1 


z 2 


1.244 


S2xpoincare 


1093 


5784 


2 


2 


2 gen. 


22.493 


S_5_100_472 (stacked sphere) 
Hom_C5_K5 (S* x S 2 ) 


480 


1430 











2.380 


24751 


143900 











5:22:22.520 


_HP2 


91 


455 











0.356 


contract ible_vertex_tiomogeneous 


1231 


12380 











46.431 



27 



Table 6: Distribution of fundamental groups for random 2-complexes with 25 vertices 
(100 runs for each p). 



V 

1 





V / 


F(2) 

\ } 


F(3) 


FU) 


F(5) 


F(6) 


F(7) 


F(8) 


F(9) 


\ J 


F(> 11) 


'non-free' 


0.40 


100 






































0.39 


100 






































0.38 


99 


1 



































0.37 


100 






































0.36 


98 


2 



































0.35 


96 


3 


1 
































0.34 


99 


1 



































0.33 


97 


3 



































0.32 


95 


5 



































0.31 


97 


3 



































0.30 


97 


3 



































0.29 


88 


10 


1 


1 





























0.28 


83 


16 


1 
































0.27 


83 


15 


2 
































0.26 


73 


19 


8 
































0.25 


73 


21 


3 


2 


1 


























0.24 


66 


27 


6 


1 





























0.23 


39 


39 


20 


1 


1 


























0.22 


31 


39 


20 


8 


2 


























n 01 

U.Zl 


OQ 
ZG 


0/1 
z4 




j 


( 


o 

z 


1 
1 


u 


1 
1 


U 


u 


n 


U 


0.20 


16 


33 


26 


16 


5 


3 


1 




















0.19 


8 


18 


24 


16 


18 


7 


7 


1 





i 











0.18 


9 


13 


24 


18 


19 


8 


3 


1 


4 


i 











0.17 


2 


6 


10 


20 


9 


21 


13 


10 


5 


3 





1 





0.16 





1 


3 


13 


12 


17 


21 


9 


6 


4 


6 


8 





0.15 











4 


16 


5 


10 


14 


10 


10 


7 


24 





0.14 

















1 


3 


5 


11 


13 


5 


62 





0.13 




















1 


1 





3 


6 


89 





0.12 
































3 


96 


1 


0.11 



































92 


8 


0.10 



































98 


2 


0.09 



































99 


1 


0.08 



































100 





0.07 



































100 





0.06 



































100 
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